Superconducting Transition in Doped Antiferromagnet 
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We investigate the superconducting transition in a doped antiferromagnet. Based on the phase 
string framework of the t-J model, an effective model describing the phase-coherence transition 
is obtained and is studied through duality transformation and renormalization group treatment . 
We show that such a topological transition is controlled by spin excitations, with the transition 
temperature determined by a characterisitic spin excitation energy. The existence of an Ising-like 
long range order of staggered current loops is also discussed. 
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One of puzzles in high- Tc cuprate superconductors is 
that the behavior of the energy gap in the quasiparticlc 
channel is quite different from that of the superconduct- 
ing transition temperature Tc- In the underdoped regime, 
two scales show divergent doping dependence: the en- 
ergy gap seems to increase as the hole concentration S 
is reduced, while Tc itself monotonically decreases [Q. 
The former is even present above Tc- By contrast, in the 
BCS theory of superconductivity, the quasiparticlc en- 
ergy gap is simply proportional to Tc at T = and van- 
ishes above Tc. So Tc seems controlled by a rather differ- 
ent low-energy physics in the cuprate superconductors. 
For example, the well-known Uemura plot Q shows the 
direct proportionality between Tc and the phase stiffness. 
Emery and Kivelson have conjectured ||], for underdoped 
cuprates, Tc is decided by the phase coherence of the 
pairing order parameter. On the other hand, a quan- 
titative feature has been recently revealed by inelastic 
neutron-scattering experiments. It was found Q that Tc 
is correlated, in roughly a linear relation, with the char- 
acteristic spin energy scale Eg associated with a mag- 
netic "resonance-like" peak ||^,^, which decreases with 
the doping concentration. 

How to establish the connection of Tc with either the 
phase coherence or the spin resonance energy [0J|] 
is currently a hot topic in various modelling, but few at- 
tempt has been made to put both the phase coherence 
and spin resonance energy within a single framework in 
understanding the mechanism of superconducting transi- 
tion in the cuprates. In this paper, we approach this issue 
by using a microscopic theory of doped antiferromagnet 
based on the t—J model. In such a description, one starts 
from the half-filling where the antiferromagnetism is well 
understood. The doping will then introduce the so-called 
phase string effect as doped holes pick up sequences of 
nontrivial signs from the spin background during their 
hopping Such a phase string effect at finite doping 
will destroy the long-range antiferromagnetic spin order, 
leading to a sharp "resonance-like" peak at a doping- 
dependent energy scale Eg ~ JS in the superconduct- 
ing phase, characterized by the holon Bose condensation 



Due to the same phase string effect, spin excita- 
tions in the superconducting state induce strong phase 
frustrations on the holon concentration and eventually 
destroy the phase coherence of the latter at a finite tem- 
perature. We will demonstrate that Tc as the onset of 
the phase coherence temperature is indeed proportional 
to Eg, which is not directly associated with the "energy 
gap" in the quasiparticlc channel. As a by-product of the 
microscopic theory, we also show that there exists stag- 
gered current loops in the superconducting ground state, 
with a hidden broken Z2 symmetry. 

We start from the t-J model. In the slave-particle 
representation one can take the aforementioned singu- 
lar phase string effect into account by a decomposition 



where h] and bia are bosonic holon 



and spinon operators, respectively. Using it, the local 
singularity of phase strings can be "gauged away" from 
the Hamiltonian and kept in the phase factor e^^'" , while 
the long distance topological properties are explicitly re- 
tained in the Hamiltonian It is believed that the re- 
sulting nontrivial topological effect of phase strings is the 
key to understanding the physics of the two-dimensional 
(2D) t-J model. Based on the exact reformulation of 
the t — J model with explicitly incorporating the phase 
string effect, an effective Hamiltonian was obtained [0: 
Hcff = Hh+Hg, in which the holon Hamiltonian is given 
by 



h\hj + H.c. 



(1) 



<tj> 



Here the phase string effect is precisely tracked by the 
lattice gauge field Afj and satisfying J^a^tj = 
'^/'^T,ienT,a'^'^ia and J2n4'zj = ±7^ per plaquette, re- 



spectively. Note that denotes the spinon number op- 
erator at site I. So Af,- will mediate the main influence 
of the spinon degrees of freedom on the holon part in a 
form of gauge field. 

In this framework, the superconducting phase is real- 
ized by the holon Bose condensation ||l^. To study its 
phase transition, we shall assume that the amplitude of 
the holon condensation has been formed at some charac- 
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teristic temperature T* > Tc such that the holon opera- 
tor can be written as hi = y^p^^e*^' {pt — 5). Then the 
holon Hamiltonian can be reexpressed as 

Hh = -PhgY. cos[0, - % - <P% - Ay , (2) 

where g = 2th. In the spinon resonating- valence-bond 
(RVB) background, the paired spinous will not directly 
contribute to the lattice gauge field A^^ according to its 
definition. So only thermally excited spinons will be seen 
by holons through Al^ in (|^). Without the lattice gauge 
field A",,, the holon I3ose condensation would occur as 
a conventional Kosterlitz-Thouless (KT) transition (0° 
as a TT flux per plaquette mainly introduces additional 
frustrations to such a system). The corresponding tran- 
sition temperature Tkt has been estimated to be 
over lOOOiir, using the parameters of the t — J Hamilto- 
nian, which is about one order of magnitude higher than 
observed in the cuprate superconductors. In the follow- 
ing, we will show how Al^ can effectively bring down the 
transition temperature to a value determined by the spin 
characteristic energy, consistent with the experiment. 
We note that the spinon part is governed by 

Hs = -JsE<^,>.e'^^'^blb]-. + H.c. At half-filling, 
without e^'^^'i , Hg reduces to the mean-field version of 
Schwinger-boson representation of the Heisenberg model 
, which well captures the antiferromagnetic correla- 
tions there. Upon doping, X^n^ij- = l/^Sien'^/' ("^f 
is the holon number operator) describes fictitious quan- 
tized TT fluxoids bound to holons and seen by spinons. 
So A^j will represent the doping infiuence on the spinon 
part, a consequence again due to the phase string effect 
1^. In the superconducting phase as holons are Bose- 
condensed, Hs has been studied in Ref. 9, where it 
is shown that a resonance-like peak will emerge in the 
spin dynamic susceptibility at (tt, tt) at a finite energy 
Eg = 2Es with Es denoting the single spinon excitation 
energy. If higher energies are neglected at low T, one 
has Hs « Eg/2J2scrlLlso- where 7so- is the Bogoliubov 
operator for spinon excitations jl^. Note that A'lj only 
depends on the density of holons. So one expects that 
the peak is still present at Tc < T < T* as long as the 
amplitude of the holon condensation persists. To first 
order approximation, in the following we shall treat Eg 
or Es as a T-independent quantity below T*. 

Now we can write down the partition function corre- 
sponding to the present system as follows: 

Z^Y. I DOD'AlD'cj^l X 

f3j2 cos[A^e(r)-0°(r)-A;(r)]-|a «2„ 

(3) 

where /? = Phg/T , lij^ = jla-Jscr^ ^^'^ primes in 
D'A^^ and D'cj)^^ imply that A^^ and (/)° satisfy the con- 



straints on X^n^^ ^^'^ Sn'/'/l given above. For conve- 
nience, the subscript p is introduced here to denote the 
link ij, and 6i — 6j is replaced by A^0(r). 

First of all, let us discuss the role of 4''^{r). Note 
that the Hamiltonian (|^) has the form of an extended 
xy model. Apart from the U{1) symmetry, there is also 
an additional local Z2 symmetry which corresponds to 
the invariance for a transformation X^a^^i^ ~ ^''^ ~^ "F"" 
at each plaquette (one can realize this by changing a 
link phase by 27r within each plaquette). Like the xy 
model, one may reexpress the partition function in 
the Coulomb gas representation through a standard du- 
ality transformation In this representation there 
are three different topological charges on the dual lattice 
site, corresponding to the vorticities of a supercurrent 
loop induced by the tt ffux of (t>'^{r), vortices bound to 
spinons through gauge field Afj, and the conventional 
2tt vortices, respectively. All these topological charges 
are coupled with each other through long range logarith- 
mic interactions. This lead to correlations among the 
topological charges. It is clear that topological charges 
with opposite vorticities have a tendency to pair at low 
temperature. Especially the pairing between the oppo- 
site charges of topological vortices related to 0° means 
that local Z2 symmetry can be broken at low tempera- 
ture jl^. One can expect that apart from the topolog- 
ical transition, there is also a Ising-like long rang order 
of staggered current loop at low temperature Be- 
cause of there are two degenerate ground states and, 
correspondingly, there are two low energy modes [ p^ . 
By constructing a Landau-Ginzburg- Wilson description, 
one can demonstrate that (j^) may be decomposed into 
two coupled extended xy models. In this representation 
the topological transition is determined by the extended 
xy terms corresponding to each low-energy mode, while 
the Z2 symmetry broken is determined by the renormal- 
ized behavior of the coupling constant h between the 
two modes where h is renormalizcd to strong coupling, 
and the two modes are locked with relative phase or 
TT |]l5| , |l6| . Since generally the Z2 symmetry broken tem- 
perature is higher than the topological transition temper- 
ature , in the following we mainly focus on the low-T 
phase without further considering the fluctuation effects 
induced by 0° (r) . The main effect of the modes coupling 
due to (t>'^{r) will be represented by the renormalization 
of /3. One can introduce an effective (3' = p^g' /T with g 
being replaced by g' to denote the effect. A further study 
on the Z2 symmetry-broken and staggered current loop 
is to be given elsewhere. 

The distinctive feature under this single mode approx- 
imation is the spinon-vortex introduced by A* with the 
corresponding vorticity given by n(r* )7r = ei/^A^A^(r) = 

7r/4X/gn J2cr ^''^\(T I^y means of the standard Villain ap- 
proximation |l8[ and duality transformation p^ ] , we can 
then arrive at the following form 
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(4) 



where <^ is the "spin wave" variable of the xy model, 
m{r*) = 0,±1,±2,..., represent the topological charges 
of the ordinary 27r vortex of the xy model, and n{r*) 
represents the topological charge of the tt vortices asso- 
ciated with spinon- vortices, where r* denotes a lattice 
site dual to r. From we see that there are two types 
of vortex perturbations for the '"spin wave" fixed line of 
the system: one is the vortex with vorticity 27r and the 
other is the vortex with vorticity tt. According to Amit 
et. al ligj] , the perturbation of vortex with vorticity wk 
has the scaling dimension n? /2 near the critical point. 
With this result we see that the 27r vortex perturbation 
is irrelevant compared to the tt vortex. It will not affect 
the critical behavior of the system. We can thus only 
focus on the effect of the tt vortices, associated with ex- 
cited spinous, in the following study of low-temperature 
topological transition. The last term on the right hand 
side of (||) controls the fluctuations of n(r*). In fact, one 
can write down the relation: 



n{r*) ~ ^|u;,(r*)|Vn2CT 



(5) 



in which Ws{r) is the single-particle wave function for 
spinous at Eg and in the coherent-state represen- 
tation |20| Ws{r) ~ Cse~''^~^<'^ /'''o, with the nor- 
malized constant. In this representation the quantum 
numbers of the states are denoted by the position of 
their centers which form a von Neumann lattice with 
lattice constants as = hs — 2lo which is a function of 
the density of doped hole: Iq = At low tem- 

perature r <C Eg, the terms with nj^ = domi- 
nates, and the terms with n^g. = 1 can be regarded 
as small corrections. Then ^ exp{7ri ^ n(r*)0(r*) — 

Ej2TY.nU ^ exp{yfcEcos[7rE I ^s(r*) ? (/.(r*)]}, 

where yt = 1/ cosh.{Eg/2T) (j/f, << 1 ), in deducing it all 
the terms with order 0{yj^) (n > 1 ) have been omitted. 
It is easy to understand that yt measures how easy (dif- 
ficult) to excite a spinon. Finally the continuum form of 
the partition function is obtained as follows (for conve- 
nience in the following we use r instead of r* to denote 
the dual lattice site): 

/ d^r\w,{r)\^-4>{r)]-y /dV,[2i^ / dVS, |«,, (r) 1^ ^(r)]^ ^q-^ 

where a is the lattice constant. We have introduced a 



tition function, which will be generated by renormal- 
ization group (RG) procedure discussed below with its 
initial value being zero. The partition function is now 
parametrized by three quantities: (3' , y, and t/6. 

Compared to the conventional xy model, the partition 
function (|^) looks more complicate. The physical ori- 
gin of this complexity comes from the spinous, which 
are at centers of the tt vortices. Note that an excited 
spinon does a cyclotron motion [ p^ , and in our coherent- 
state representation its distribution function | Wsir) p 
has an attenuation radius Zq- So the tt spinon- vortex 
has a finite vortex-core of radius Iq, which is reflected 
in (H) through | Ws(r) p. We shall treat the prob- 
lem by means of the Wilson's RG analysis |^ , p2| . 
In this RG procedure we first divide (/) into the high 
energy part and low energy part, respectively: (/)a = 
+ h, where 0A(r) = J^^^_^j^d'^p/ {2Tr)^ (j){r) exp[ip.y:], 
0A'(r) = /o<p<A' dV(27r)20(r)exp[zp.x], A' A - dA, 
and A = 1/a is the momentum cut off; then the high 
energy part h is averaged out by means of the cumu- 
lant expansion; finally we rescale the system so as to 
restore the original cutoff. Because of the presence of 
a characteristic length scale 21q of the vortex core, the 
RG analysis is separated into two steps; we first treat 
local physics within the length scale 2lo, then study the 
low energy and long wave-length physics beyond such a 
scale. After some algebra, following recursion relations 
can be obtained 



dlo — —Iq — , 
a 

P' 1 



dy. 



da 



4 (AZo)2 
dy = ci(3'^yl 



Vb- 



2 2 1 da 



[Alo] 



(7) 
(8) 
(9) 



where ci = O.OIStt^. From (0)-(||), the effective param- 
eters at the length scale 2lo can be determine: y^ = 

y(A;) = 2cif3' \^l^..^^^^^^ , 



2/fc(A/) 



V 



cosh(Es/T) 

where A/ = l/(2Zo)- With the length scale 21q being 
reduced to our "new" lattice constant a, the effective 
attenuation radius of | Ws{r) p also "shrinks" from Iq 
to a/2. So the effective range of | Ws{r) p is within 
the unit plaquette of the new lattice now, which can 
be reasonably treated as the S{r) function, such that 
J d^r I Ws(r) p (/)(r) ~ 0(fs)- The kinetic term in 
can be then rewritten as — 1(1 + 2y') / d'^r[d(l>{r)Y . 



After the variable change i 



6, we finally obtain 



term y J d^r^^ / d'^rds \ Ws{r) 



into the par- 



where /3" = f3'/{l + 2y'). Equation (|l^) is exactly the 
partition function of sine-Gordon model [ p^p^ , |2^ with 
effective parameters and f3" . 
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(c) 
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FIG. 1. (a) Tc as a function of the spin resonance-like en- 
ergy Eg at different dopings, from top to bottom correspond- 
ing to 5 =0.2, 0.15, 0.1, and 0.5 respectively, (b) Tc versus 5 
(solid) with using Eg{S) shown in the inset and g' = 0.2eV. 
The dashed line represents the ordinary KT transition tem- 
perature without considering spinon- vortices, (c) Eg{5) ver- 
sus Tc{S) is ploted for the same parameters as in (b). The 
Inset: Eg{5) versus 5 obtained in Ref. p4| . 

This model has been well studied, and the transition 
temperature (i.e., Tc), at which pairing between spinon- 
vortices and -antivortices dissolves, is determined by the 
following equation (§21: 



tively brings Tc down to the right order of magnitude 
as the consequence that the spinon-vortices instead of 
the conventional 2tt vortices control the superconducting 
phase coherence transition. In Fig. 1(c), Eg{5)— Tc{5) 
is shown, which also are both qualitatively and quanti- 
tatively in good agreement with the experimental results 

i- . 

In conclusion, we have established a quantitative the- 
ory of superconducting transition based on an effective 
spin-charge separation description of the doped AF Mott 
insulator. The underlying physics is that the phase coher- 
ence transition is controlled by thermal spin excitations, 
which substantially reduce Tc from Tkt to a fraction of 
Eg/kB. It resolves the long-standing issue of how Tc can 
be quantitatively connected to the characteristic spin en- 
ergy scale in a doped AF Mott insulator. The obtained 
Tc — S and Eg — Tc relations are in good agreement with 
those observed in the cuprate superconductors, lending 
a strong support for the experimental relevance of the 
spin-charge separation theory. 
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2 = yU2. 



(11) 



Numerical results are shown in Fig. 1. Tc as a, function 
of Eg is plotted at several doping concentrations in Fig. 
1(a). It shows that Tc increases monotonically with Eg 
and saturates at larger Eg. According to Fig. 1(a), Tc/g' 
linearly scales with Eg/g' at small Eg/g' where Tc/Eg 
does not depend on g', with the ratio (^ 1/4 — 1/5) 
being only weakly doping dependent. So Tc essentially 
is determined by the characteristic spin resonance-like 
energy Eg. For instance, plugging in the experimental 
value of Eg = 41 meV for the optimal doping YBCO 
compound one estimates Tc ^ 100 K, very close to 
the experimental value. Since Eg as a function of S has 
already been obtained in the same framework ||2^ as re- 
plotted in the inset of Figs. 1(b) and (c) (with Jg = 0.1 
eV in Ref. |Q), one can use such a calculated Eg{5) to 
determine Tc vs. S as shown in Fig. 1(b) (solid curve). 
For comparison, the dashed line in Fig. 1(b) represents 
Tkt (— 7r(5(7'/2) without including the spinon-vortices 
due to A^^. Here we choose g' ~ 2th — 0.2 eV, and find 
T, 



KT 



y . Here we choose g' ~ zt/i 
510 K at S = O.U while n = 107^. As noted 



before, Tc is not sensitive to g', but Tkt does. So if 
we take g' = 0.5 eV at the same S = 0.14, Tc increases 
to UAK but Tkt reaches 1276 K. Therefore, Eg effec- 
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